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This article, the first of two, describes and analyzes the thermal performance of a 
temperature control system for the antenna microwave transmitter (klystron tube). In 
this article, only the mathematical model is presented along with the details of a com- 
puter program which is written for the system simulation and the performance param- 
eterization. Analytical expressions are presented in Appendixes. 


I. Introduction 

The microwave transmitter is one of the vital subsystems in 
the Deep Space Network Antenna Tracking System. The trans- 
mitter (klystron tube) converts the electrical energy to a 
microwave signal with a conversion efficiency in the order of 
40-60 percent. Due to the inherently large heat generation in 
the transmitter and the sensitive temperature control require- 
ments, a well-designed cooling system is required. The trans- 
mitter cooling design is commonly made, based on past engi- 
neering experience taken from the working transmitters. A 
detailed thermal analysis of an optimized system was not 
necessary. However, due to the increasing complexity of add- 
on components in the electrical circuit, for better controls and 
continuous upgrade of the existing system performance, the 
need for a new cooling system with an optimum design is 
found inevitable. In order to support this new cooling system 
design, this work is initiated to give a detailed thermal analysis 
and to evaluate the thermal performance of the system. A 
short computer program was written to analyze the transient 
cooling process. This article, the first of two, includes the 
system governing equations needed for a detailed transient 
thermal analysis and a detailed description of the simulation 
program. The second phase analysis is intended to evaluate the 


thermal behavior of the new system under varying operating 
conditions, with different component selection, location and 
size on the temperature control requirements. 

II. System Description 

Figure 1 illustrates one proposed design of the new micro- 
wave transmitter cooling system. The numbers assigned to the 
various points in Fig. 1 represent the fluid stations. The system 
is comprised of (1) a cross flow air-to-liquid heat exchanger 
(A), (2) a temperature control valve (B) which controls the 
amount of fluid flowing out from the heat exchanger (sta- 
tion 1) and the bypass of the returning flow (station 3) to 
deliver the desired fluid temperature (station 4), (3) a filter (C) 
which removes contaminants resulted from the corrosion in 
the system; and (4) a pressurized storage tank (D) which per- 
forms as a thermal flywheel or damper for all possible thermal- 
fluctuations in the circuit. The storage tank is charged with 
nitrogen gas for maintaining the desired pressure. An electric 
resistance heater is provided inside the storage tank which 
operates only when the control valve (B) fails to supply a fluid 
temperature above the temperature setpoint of the tank. A 
pump (E) pumps the fluid from state 9 to 10 and maintains 
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the fluid in the loop at the design pressure. A flow meter (F) 
monitors the system flow rate at station 1 1 . A purity loop (G) 
purifies the working fluid of the cooling system from state 31 
to state 32. A pressure control valve (H) controls the working 
pressure at the microwave transmitter station 17. A flow 
bypass valve (I) regulates the proper flow to the microwave 
transmitter. The microwave transmitter includes two major 
assemblies, a water load (J) and the klystron assembly itself 
(K). The changes in fluid states due to the piping pressure drop 
and heat losses are given different numbers as in the stations 
1-2, 4-5, 6-7, 8-9, 10-11, 12-13, 13-14, 13-31, 15-16, 16-17, 
16-26, 17-18, 17-21, 19-20, 22-23, 27-28, 32-33,34-35,35-3, 
and 35-36. 


III. Thermal Analysis 

The following assumptions were made in the mathematical 
formulation of the system: 

(1) The entire system, excluding the storage tank, is 
assumed at steady state, with the fluid flowing at a 
constant rate. 

(2) The system is located in an environment with uniform 
ambient temperature and solar radiation. 

(3) Axial conduction heat transfer is neglected from one 
end of a component to another. 

(4) Sky and ambient temperatures are assumed approxi- 
mately the same for simplicity. 

(5) The fluid inside the storage tank is assumed to be 
thoroughly mixed to give a uniform temperature. The 
fluid temperature leaving the tank is assumed to be the 
same as the fluid temperature inside the storage tank. 

(6) The walls of the tank and the various pipes are assumed 
isothermal. 

(7) The temperature of nitrogen gas inside the storage tank 
is assumed to be the same as the fluid temperature, 
with negligible heat transfer effect. 

(8) The storage tank is assumed in the form of a cylinder 
exposed partially to solar radiation. For simulation 
purposes, the surface area of the storage tank which is 
exposed to solar radiation is assumed to be a ratio X of 
the total tank surface area. 

Using the above assumptions, two systems, an insulated 
system and a noninsulated system, have been analyzed. Appen- 
dixes A and B give the details of the heat balance for the 
insulated and noninsulated systems. Tire temperature distribu- 
tion of any piece of piping, whether it is insulated or not, is 
characterized by the expression in Eq. (A-25). 


T f (x) = T fe + 


' E 1 - *1 (Tf, e - T a)\ (] 
B l _T 


~c,x 

e 1 ) (1) 


where B j is the overall heat transfer coefficient between the 
fluid and the ambient air, Ey is the solar energy absorbed by 
the pipe, Tf e is the fluid temperature to the pipe entrance, Cy 
is the pipe characteristic constant, x is the distance measured 
from the pipe entrance, and 7^(x) is the fluid temperature at 
that location. 


The general solution for the time history of the storage 
tank outlet temperature can be expressed as in Appendix B by 
an equation similar to Eq. (1) as, 


T f (t) = 7^(0) + 


D 2 


-D.t 
(1-e 1 ) 


( 2 ) 


where 7y(f) is the fluid temperature leaving the tank after a 
time internal t, and D y and D 2 are the storage tank character- 
istic constants that depend on the mass flow rate, thermal 
conductances, and the tank inlet fluid temperature as shown in 
Appendix B. The value D 2 /Dy represents the final equilibrium 
temperature of the tank with a constant inlet temperature. 

In order to determine the time increment for a stable and 
convergent solution of the transient analysis, the time constant 
1 /Dy is calculated. When time t equals Dy, the temperature 
difference (D 2 lDy)~ 7y(0) would drop by 36.8 percent from 
its initial value. A new time increment is calculated for each 
time interval. 

The temperature rise due to the pump dissipating work can 
be determined by using the thermodynamics relationship for 
the incompressible flow as 


W_ 


AT f = m f C p 


( 3 ) 


where W p is the actual pump work in the system and my- and 
C p are the fluid mass flow rate and the fluid specific heat, 
respectively. The heat dissipation in the system can be corre- 
lated with the temperature rise of the cooling fluid. The 
general expression for the fluid outlet temperature T^ ex is 


T f,ex 


Q + m fCp T f,e 

™f C P 


( 4 ) 


where Q is the heat generated by the klystron systems and Tj- e 
is the fluid temperature at the klystron inlet. 
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Stations 4, 24, 29, and 34 represent the fluid mixing 
locations. The general expression for the mixed fluid tempera- 
ture is 


T r,m + 


( 5 ) 


where P is the mass ratio of the flow rate in the first stream, 
whose temperature is 7^ t , to the total mass flow. Tf 2 is the 
fluid temperature of the second stream. The fluid temperature 
at station 35 represents the return fluid temperature just 
before the heat exchanger bypass loop. The bypass ratio y is 
determined by the temperature control valve and is calculated 
by iteration at each time step. The final y value is determined 
when the system reaches the steady state operating conditions. 
Due to the insufficient physical information of the heat 
exchanger, it is assumed that the heat exchanger is of the 
counter flow type. The fluid temperature leaving the 
exchanger, Tj- ex , can be expressed in terms of the entering 
fluid temperature Tf e as (Ref. 1) 


1 f,ex 


T - 
f,e 


a>, e - yv-en 


K C a 


- P P 


( 6 ) 


where 


P 


UA 




UA = the product of the overall heat transfer coefficient 
and the total surface area given in Ref. 2. 

C a = the specific heat of the air 
C p = the specific heat of the fluid 

= the mass flow rate of the air 
riij- = the mass flow rate of the fluid 


IV. Computer Program Methodology 

A computer program entitled Kystron Cooling Project 
(KCP) is written using the mathematical model developed in 
Section III and Appendixes A and B. The program includes a 
main program and five subroutines. The program starts its 
computation at the heat exchanger outlet (station 1). Initially, 
the fluid temperature is assumed 10°F above the ambient 
condition and all fluid is flowing through the heat exchanger 
(7 equals to zero). The program then proceeds its computation 
in the sequence as described in Fig. 1. For each of the pipe loss 
calculations, the subroutine PIPE is called to compute the 
fluid exit temperature. The storage tank calculation is handled 
in the TANK subroutine, and the time constant (1/Dj) will be 
used for the next time increment. In each of the pipe loss 
calculations, HTCOE and AHTOE subroutines are called to 
calculate the heat transfer coefficient of the fluid inside the 
pipe and the heat transfer coefficient to the ambient air, 
respectively, for different pipe sizes and lengths. The 
WRTEMP subroutine prints out the fluid condition at each of 
the stations indicated in Fig. 1 . A temperature sensor is built 
into the program which compares the calculated temperature 
at station 18 and the design setpoint of the klystron assembly. 
The difference of the two temperatures will be adjusted by the 
temperature control valve at station 2. The valve determines 
the amount of the fluid flow passing through the heat 
exchanger which mixes with the returning flow (3) to deliver 
the new adjusted temperature by iterations. When the design 
condition of the klystron assembly is met for three consecu- 
tive times, the system has reached steady state operating 
condition, and the time required for the first occurrence is 
considered to be the system response time. 


V. Summary 

The computer program has not only analyzed the perfor- 
mance of a cooling system, but it also has provided a helpful 
tool for the designer in setting design requirements. However, 
in order to evaluate the performance of a system, detailed 
information of the system is needed. When this study was 
initiated, several cooling design concepts were considered. 
Physical dimensions and system control requirements were not 
available. The second phase of this study will report the 
performance and design requirements of a finally selected 
cooling system. 
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Definition of Terms 


A area 

B l - /? 6 thermal conduction 

C a air specific heat 

C p fluid specific heat 

- C 2 constants 

D diameter 

Z)j - £>2 constants 

Ey - E 4 energy flux 

G heat capacity 

Hy - H 1 heat rate 

h convection heat transfer coefficient 

/ solar flux 

J constant 

K thermal conductivity 

L pipe length 

M total fluid mass inside the storage tank 

m mass flow rate 

Q heat rate 

T temperature 

t time 

U effective heat transfer coefficient 

V velocity 

x distance 


Y thickness 

a absorptivity 

p density 

p viscosity 

5 parameter 

e emissivity 

(3 mass flow ratio 

7 heat exchanger bypass ratio 

X ratio of the tank surface area exposed to insulation 

7, surface area ratio of tank wall and tank insulation 


Subscripts 


a 

ambient air 

e 

entrance 

ex 

exit 

f 

fluid 

I 

insulation 

i 

inside 

m 

mixing fluid 

o 

outside 

ti 

tank insulation 

tw 

tank wall 

w 

pipe wall 

X 

location 
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Appendix A 

Derivation of Temperature Distribution for a Pipe 


In this appendix, the derivations of the heat transfer equa- 
tions for both an insulated and a noninsulated pipe are made. 
Following the assumptions made in Section III, Fig. A-l illus- 
trates a segment of an insulated pipe whose length is dx 
located at a distance x from the inlet fluid section. The 
differential rates of heat flux are divided as follows: 

dQ x = total absorbed solar insolation (direct and diffuse) 
on the outer insulation surface 

dQ x = otjl D J o dx (A-l) 

where aj is the absorptivity of tire pipe insulation 

dQ 2 = effective heat transfer between the insulation outer 
surface and the ambient air which combines the 
convection and radiation parts 


Writing the energy balance equations for each of the elemen- 
tary pipe components at steady state will yield the following 
equations: 

For the insulation surface: 

dQ x - dQ 2 -dQ 3 =0 (A-6) 

For the pipe wall: 

dQ 3 - dQ 4 = 0 (A-7) 

For the fluid: 

dQ 4 ~ dQ 5 = 0 (A-8) 


d Q 2 = U l0\ o -T a )nD [o dx (A-2) 


where U { is the effective heat transfer coefficient 

dQ 3 = conduction heat transfer from the insulation layer 
to the pipe wall 


dQ, 


C T I,o- T W ? dx 
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2 ttK.. 


(A-3) 


dQ 4 = convection heat transfer between the interior pipe 
wall and the fluid flowing through with a heat 
transfer coefficient hf 


The coefficient Uj given in Eq. (A-2), defined as the effec- 
tive heat transfer coefficient between the insulation outer 
surface and the ambient air, can be expressed as the summa- 
tion of a convection heat transfer coefficient h a and a “lin- 
earized” radiation heat transfer one: 

U I = h a + 6e (7/, 0 + ^a) (7/,o + T'a) (A-9) 

The coefficient h a can be determined by the expression 
(Ref. 3) 


h 


a 


am 


(A- 10) 


where J and n are constants depending upon the Reynolds 
number as given in Ref. 3 and V a is the air velocity. 


= T f* nD wj dx (A- 4 ) 

dQ 5 = sensible heat carried by the fluid from section x to 
section (x+c/x) 


The coefficient hf given in Eq. (A4), defined as the convec- 
tion heat transfer coefficient between the pipe inner surface 
and the fluid, is obtained from the Nusselt number Nu D and 
the diameter D w ,• as follows: 


dQ s = m f C p 



(A-5) 


Nu n K, 
h = — — — — 
f D . 

w,i 


(A-l 1) 
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The Nusselt number can be determined by applying any one of 
the following expressions depending on the flow condition: 

(a) If Reynolds number Re D is greater than 7000, the 
flow is fully developed turbulent (Ref. 3) and the 
Nusselt number is given by 


Nu d = .023 Re D - 8 Pr 3 (A-12) 

where Pr is the Prandtl number 

(b) if Re is less than 2100, the flow is laminar and the 
Nusselt number is obtained from 


Nu d = 3.66 + 


Ml^f) R e D 


Pr 


1 + 


- 04 [(%)*'* 


12/3 


Pr 


(A-13) 


where L is the length of the pipe 

If Reynolds number ranges 2100 < Re D < 7000, the flow 
is considered as a transition flow. Since very few experimental 
data are available in this flow regime, a straight line inter- 
polation for the Nusselt numbers between the laminar regime 
and the turbulent regime is assumed as a first approximation 


Nu d = Nu d (2100) + (Re D - 2100) 


'Nu d (7000)- Nu d (2100)“ 


7000-2100 


(A- 14) 


Equations (A-6) through Eq. (A-8) can be rewritten after 
substituting all the dQ values and dividing by (D/ 0 dx) as 
follows: 

For the insulation surface 


E 1 - B X <- T I,o - V “ *2 c Tj,o - T W)f ) = o (A-15) 
For the pipe wall: 

V^o-W-M^.r 7 }^ 0 < A - 16 ) 


For the fluid: 



Expressing the temperature 7} 0 and T w ,-in terms of T a using 
Eqs. (A-15) and (A-16) yields 


E l+Bl T a+ B 2 T wJ 

B a 


E 2 + B 3 T f 


w,i B c 



where 
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Substituting T w ,• in Eq. (A-17), the differential equation 
for the fluid in the pipe at any position x from the entrance 
will be 


dT f 

dx 


+ C\Tf 


(A-22) 


The quantity [ZTj - By {T^ e - T a )] represents the net energy 
exchange between the fluid and the ambient air. By is, there- 
fore, the overall heat transfer coefficient as given by Eq. 
(A-18), and iTj is the solar energy absorbed by the system as 
given in Eq. (A-19). 


where 



The second term of Eq. (A-25) (1 - e -c,JC ) is considered a 
dimensionless flow factor. Cj as given by Eq. (A-23) is a 
function of the heat transfer coefficients ( B’s ) and the mass 
flow rate (my). 

(A-23) 

Equation (A-25) can be applied to a noninsulated pipe, as 
well with the following changes: 


Note that Cy is positive since ( B s - B 3 ) can be proved to be (i) The absorptivity ot w of the pipe wall should be used to 

equal to ByB^jB/y. replace of the insulation in Eq. (A-l). 


The general solution of the differential Eq. (A-22) is 


C.x C 2 
T f (x) = Se 1 +^r 


(A-24) 


where 5 is an arbitrary constant which is determined by the 
following boundary condition: 

At the pipe inlet section ( X = 0) the fluid temperature 7y e 
is assumed given. Substituting in Eq. (A-24), the temperature 
distribution Tf{x) can be reduced to: 

pi - fit (T f - T )1 - c.x 

W T r.'*l - V* ') 

(A-25) 


(2) The outside diameter D w 0 of the pipe should replace 
Dj a of the insulation. 

(3) The outer surface temperature T w 0 of the pipe should 
replace Tj 0 of the insulation. 

(4) The heat conduction term B 2 in Eq. (A-18) should be 
modified as follows: 


B , 


“ 

ID VI 

£n 

n w -°) 

v w,i / 


J 


(Pw.o) 
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Appendix B 

Derivation of Temperature Variation for a Storage Tank 


In this appendix, the heat transfer equations for an dTj- 

insulated and a noninsulated fluid tank are derived, using the #5 = M (B-5) 

assumptions listed in Section III in the text. 

where M is the total fluid mass inside the tank 

/fg = extracted sensible heat gain by the fluid passing 
through the storage tank. 


Hy = total absorbed solar radiation (direct and diffuse) on 
an area with a ratio X of the outer tank surface area. 


H 6 =’ h f C pVf- T f,e) ( B ' 6 ) 


Figure B-l illustrates a cross section of an insulated storage 
tank with an insulation thickness Y). The heat flux is divided 
as follows: 


H i = ( B -D 

where a tj is the absorptivity of the insulation enve- 
lope. 

H 2 = effective heat transfer between the insulation outer 
surface and the ambient air. 

»i * w-nKi,. < B - 2 > 

where U ti is the effective heat transfer coefficient 
which is determined by equations similar to Eqs. 
(A-9) and (A- 10) and T ti 0 is the outside tank insula- 
tion surface temperature. 

// 3 = conduction heat transfer through the insulation layer 
with thickness Y tj . 



H 4 = convection heat transfer from the storage wall to the 
fluid inside the storage. 

H A h f,T^ T tw,i~ T fi A tw,i O 3 - 4 ) 

where hj j is the convection heat transfer coefficient 
between the tank wall and the fluid inside the tank. 

= sensible heat gain by the fluid inside the storage 
tank. 


where rhj- is the steady mass flow rate in and out of 
the tank, and e is the fluid temperature at the 
storage tank entrance. 7\ is assumed to be constant 
in each time step calculation. 

= heat generated by the auxiliary heater located inside 
the tank. Three energy balance equations are found 
by grouping Eqs. (B-l) to (B-6) and H s as follows: 

For the insulation surface 

H x -H 2 -H 3 = 0 (B-7) 

For the tank wall 

// 3 - H A = 0 (B-8) 

For the fluid inside the tank 

H 5 =H a +H € +H 1 (B-9) 

After substituting all the values into Eqs. (B-7) to (B-9) and 

dividing by (A ( . Q ), the energy balance equations become 

For the insulation surface 

E l- B l( T ti,o-V- B 2( T ti ,o- W ( B - 10 > 

For the tank wall 
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For the fluid inside the tank 


substituting T^ v . into Eq. (B-12) yields 


where 


dT f 


mc p* =b 3 <Tt„.r T } +G Vr T f.J 


^ ti,o 


B i = U ti 


B 2 = 


y\ 


M VW 


(B-12) where 


(B-13) 


dT f 

-^r =D i T f +D 2 


b 6 b 3 




+ G 


B c 


D 


1 M C 


(B-17) 


(B-18) 


b 3 h fjy i 


?1 


tw,i 
^ ti,o 


E l X a ti 1 


ni f C„ 
r _ / P 

G ~ A — 
A ti,o 


(B-14) 


* 

4 


D , 




^5 ~ G T f,e + E A 


MC 


The final general solution for this differential equation is 


-D,t D 2 

T f (t) = 5 e 1 (B-19) 


Using Eqs. (B-10) and (B-ll), temperatures T ti and T tw .■ , . . ... . . , . 

, , lw ’’ where 5 is an arbitrary constant determined by the initial 

can be expressed as 


condition: 


T . 
tt,o 


tw,i 


E l +B 1 T a + B 2 T tw,i 

B A 


E 2 + E 3 T f 
B c 


where 84= By + B 2 


B 5 =B 2 l 1 -J) + ^3 


■©-.•ft 


b 7 t 

Z a 


At time equals zero, the fluid temperature inside the tank 
7^(0) is given. Substituting in Eq. (B-19) and rearranging 
(B-15) terms; the temperature variation Tf{t ) is given by 


T f (t) = 7/0) + 




-D.t 

(1-e 1 ) (B-20) 


The quantity [(D 2 /D 1 ) - 7/0)] can be expressed as the ampli- 
tude of disturbance which depends on the thermal conduc- 
tances ( B’s ), the heat fluxes (TTs), the flow capacitance ( G ) 
and the storage tank inlet temperature (7/). If the mass flow 
(B-16) (/ hj) is equal to zero, the quantity would reduce to become 

the static disturbant amplitude in the stagnation condition. 
The factor (1 - e~ D 1 ! ) is defined as the time decay factor 
which describes the time history of the storage tank. 
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(4) i?2 Eq. (B-13) should be modified as follows: 


Equation (B-20) can be applied to a noninsulated tank with 
the following changes: 

(1) The outside tank area (A^ Q ) should be used to 
replace the outside tank insulation area 0 ) in all 
appropriate equations. 

(2) The tank wall absorptivity (a [w ) should replace the 
tank insulations absorptivity (a n ) in Eq. (B-l). 

(3) The outer tank wall temperature (T tw Q ) should 
replace the outer tank insulation’s temperature (T ti 0 ) 
in Eqs. (B-2) and (B-3). 



(5) The effective heat transfer coefficient U should be 
recalculated using the outside diameter of the tank 
instead of the outside diameter of the insulation. 
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Fig. B-1. Cross section of a storage tank 
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